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Two decades ago, a phenomenon resembling Landau damping was described in the synchronization of glob¬ 
ally coupled oscillators: the evidence of a regime where the order parameter decays when linear theory predicts 
neutral stability for the incoherent state. We here show that such an effect is far more generic, as soon as phase 
oscillators couple to their mean field according to their natural frequencies, being then grouped into two distinct 
populations of conformists and contrarians. We report the analytical solution of this latter situation, which al¬ 
lows determining the critical coupling strength and the stability of the incoherent state, together with extensive 
numerical simulations that fully support all theoretical predictions. The relevance of our results is discussed in 
relationship to collective phenomena occurring in polarized social systems. 

PACS numbers: 05.45.Xt, 68.18.Jk, 89.75.-k 


In the forties, the Soviet physicist Lev Davidovich Landau 
(the 1962 Nobel Laureate for his theory on superfluidity) pre¬ 
dicted the damping (the exponential decrease as a function of 
time) of electrostatic charge waves in a collision-less plasma 
lillJ]. After almost two decades of controversy, the ’’Lan¬ 
dau damping” (LD) was eventually verified experimentally 
0 , and it was even argued that similar phenomena could take 
place in galactic dynamics 0]. Mathematically, LD is entirely 
due to the occurrence of fake eigenvalues caused by analytic 
continuation, which lead to an exponential decay of the elec¬ 
tric field even when the density perturbation does not. 

More recently, Strogatz et al. investigated the synchroniza¬ 
tion transition in ensembles of coupled oscillators 0, and 
found a regime (below the synchronization threshold, and for 
which linear theory foretells neutral stability) where the re¬ 
laxation to the incoherent state is indeed exponential, with a 
decaying mechanism remarkably similar to LD. 

In this Letter, we report on theoretical analysis and numer¬ 
ical simulations that demonstrate how the presence of LD is, 
actually, far more general in the synchronization route of glob¬ 
ally coupled oscillators. Without lack of generality, we start 
by considering a frequency-weighted Kuramoto Q0] model of 
N phase oscillators, in which units are coupled to the mean 
field according to their natural frequencies: 

N 

=UJ j + -jf-Yl sin ( 0 " - N - (1 ) 

n= 1 

Here 9j {utj) is the instantaneous phase (the natural frequency) 
of the jth oscillator, dot denotes a temporal derivative, and 
n > 0 is a global coupling strength parameter. The set { utj } 
of natural frequencies is drawn from a given frequency dis¬ 
tribution (FD) g(ut) with, in general, both a positive and a 
negative domain. Eq. 0 belongs to the class of the so-called 


* Authors contributed equally to this work, 
tCorresponding author: stefano.boccaletti@gmail.com 
^Corresponding author: guanshuguang@hotmail.com 


generalized Kuramoto models BUI- At variance with Refs. 
Imlori . the effective coupling nutj can be here either positive 
or negative (depending on the sign of the natural frequency 
utj), reflecting the fact that the interactions among individu¬ 
als can inherently be repulsive in real systems 1114141711 . For 
instance, it is well known that both excitatory and inhibitory 
couplings characterize the interaction structure of neural en¬ 
sembles Oil Q0, as well as similar kinds of coupling can also 
be found in social interactions. 

Eq. ([]} is fully analytically solvable, and we here focus on 
unveiling the details of several novel phenomena of LD, which 
characterize the transition to synchronization in such systems. 
In its mean-field form, Eq. (|T]| can be written as 

6j = utj + KLUjr sin(<^ — 8j), j = 1,..., N, (2) 

where r and (ft are order parameters defined by rtN : = 
jy Y^jLi ■ Flere oscillators can be, in general, grouped 
into two populations, according to their effective coupling: 
those with positive utj will behave like conformists attempt¬ 
ing to follow the global rhythm of the system, whereas those 
with negative utj will tend to act as contrarians , i.e. they will 
always try opposing the system's global beat lfl5l4l7tl . Let us 
then analyze the synchronization transition in Eq. (|T|i in the 
presence of a Lorentzian FD g(ui) = / it) / {{ut — A ) 2 + 7 2 ], 
where A and 7 are the central frequency and the width at half 
maximum, respectively, with A actually controlling the pro¬ 
portion of conformists to contrarians in the ensemble. 

In the thermodynamic limit, i.e., N -A 00 , a density func¬ 
tion p{9 , u), t ) can be defined, where p dO denotes the fraction 
of oscillators with frequency uj whose phases have values be¬ 
tween 9 and 9 + d9 at time t. p satisfies the normalization 
condition p{9, ut, t)d9 = 1 for all uj, and all t. The evolu¬ 
tion of p is governed by the continuity equation 


dp d{pv) 
dt d9 


(3) 


where the velocity v is given by v = ut + Kutrsm{(jt — 9). 
On its turn, the order parameter can be expressed as re u ‘‘ = 
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FIG. 1: (Color online) The spectra for the characteristic equation ©. For a Lorentzian frequency distribution, the continuous spectrum 
is the whole imaginary axis. The solid squares denote, instead, the discrete eigenvalues. The purple ovals in (a), (c), and (d) mark the ghost 
(fake) eigenvalues predicted by the linear theory. As discussed in the text, such ghost eigenvalues remarkably control the decaying rate of order 
parameter r(t) in the neutrally stable regimes, exactly as in the Landau damping context. It should be pointed out that the ghost eigenvalue in 
(d) could also be above the real axis, depending on parameters. For better visualization, we do not plot the other irrelevant (fake) eigenvalues. 


r e lS p{9, w, t)g(u>)du}dO, and therefore Eq. © can be 
rewritten as 

do d f 27T I + °° 

p(0' V ,t)g(ui')dw'dO']}. (4) 

For the incoherent state, po( 8 ,ui,t) = l/(2tv). A small per¬ 
turbation from that state, i.e. 

p( 6 ,w,t) = ^- + e[c(uj,t)e ze + c*(uj,t)e~ lB + t? ± ( 8 ,tj,i)\ 

(5) 

can be considered, where e <C 1 , and 1 )^( 8 ,u>,t) represents 
the higher Fourier harmonics. Substituting Eq. © into Eq. 
0 , one gets the linearized characteristic equation 

dc(u,t) . ' ^ , KtV f + °° , , , 

——— = -uv ■ c(w, t) + — J c(w , t)g(ui Jaw . ( 6 ) 

Eq. © has both continuous and discrete spectra (details 
form part of the material contained in the SI), and one can ex¬ 
tract the critical coupling strength k c as well as determine the 
stability properties of the incoherent state (which are schemat¬ 
ically reported in Fig. ©. Specifically: 

(Case 1) A > 0: the conformists prevail over the contrar¬ 
ians. k'[ = 27 /A and the spectra are schematically plotted 
in Fig. |TJa) and (b). When k > kJ, Eq. (© has a continu¬ 
ous spectrum on the imaginary axis and a discrete eigenvalue 
in the right half complex plane. Accordingly, the incoherent 
state is unstable. When 0 < k < /tj, there is no discrete 
eigenvalue, and the incoherent state is neutrally stable. 

(Case 2) A = 0: conformists and contrarians are equal in 
number. In this case, = 00 , implying that synchronization 
can never be achieved. As shown in Fig. ©c), here no discrete 
eigenvalues exist outside the imaginary axis, and therefore the 
incoherent solution is always (i.e. for any arbitrary coupling 
strength) neutrally stable. 

(Case 3) A < 0: the contrarians prevail over the con¬ 
formists. Kg = — 27 /A, and the spectra are shown in Fig. 


©d) and (e). When k > Kg there is (besides the continuous 
spectrum on the imaginary axis) an eigenvalue in the left half 
complex plane. Thus, the incoherent state is linearly stable. 
When 0 < k < Kg, there is no discrete eigenvalue, and once 
again the incoherent state is neutrally stable. 

The current situation shares connections and differences 
with the classical Kuramoto model 01 - For A > 0, the 
stability of the incoherent state is the same as that of Refs. 
0! , though the equations of the two models are essentially 
different. For A < 0, the incoherent state changes from neu¬ 
trally stable to stable when k exceeds Kg, and the bifurcation 
here only involves a change in stability of the state, while no 
new steady states emerge. Finally, for A = 0, the incoher¬ 
ent state is always neutrally stable, regardless of k. The two 
latter phenomena are novel, and inherent in our frequency- 
weighted Kuramoto model, which allows the two populations 
of conformist and contrarian oscillators to coexist. 

In order to unveil LD in the model, we analytically extract 
the equation ruling the relaxation behavior of r(t) in all three 
neutrally stable regimes predicted by the linear theory, i.e., 
the cases corresponding to Figs, ©a), (c), and (d). As it will 
appear momentarily, the most remarkable finding is that r(t) 
decays, indeed, exponentially in all cases. The use of Eq. © 
into the definition of r leads to r(t) = 27re|R(f)|, where 

/ +00 

c{t,uj)g(uj)duj. (7) 

-OO 

Then, Eq. © becomes 

dc(u ), t ) . . KW , 

— — — =-iwc(w,t) + —R{t). (8) 

For any given initial condition co(w) = c(0,w), the solution 
of Eq. ® is: 

r 4 1 

c(t,w) = e~ lult [ / -Kwi?(r)e IUJT dr + c 0 (w)]. (9) 
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Substitution of Eq. (0 into Eq. 0 leads to 

/ +oo ft 1 

e~ lult [ / -KU>R(T)e lulT dT]g{u>)du} 

-oo J o 2 

/ +oo 

e~ lult co(u])g(uj)duj 

-OO 

_ 1 /■* „ 

= (c 0 g)(t) + -Kj G(r)R(t — r)dr, (10) 

where the hat denotes the Fourier transform, i.e., (cog)(t) = 
e~ iult c 0 {uj)g(u])duj and G(t) = f™ oo e~ iult ujg(uj)duj. 
The solution of Eq. (fTOt is obtained by a Laplace transform 
(all details are reported in the SI), and reads as 

R(t) = ^ = exp{[y (« - J) - i{^ + A)]i}. (11) 



FIG. 2: (Color online) Landau damping effects in the order pa¬ 
rameter R(t). Horizontal lines is time, solid lines refer to the nu¬ 
merical solutions of Eq. 0 . and dashed lines are the analytical pre¬ 
dictions of Eq. CD. The three log-linear plots refer to the cases 
A > 0 (a), A = 0 (b), and A < 0 (c). In the numerical simulations, 
the initial states of the system are set in the fully coherent states. 


Notice that, in the derivation of Eq. CD, we do not yet im¬ 
posed specific constraints to the functional form of p,, there¬ 
fore [based on the closed form of R(t)], we here below sum¬ 
marize what happens to the relaxation behavior of r(t) in the 
three different cases of A, while addressing at the same time 
the reader to the SI material for all technical details. In the 
following, Ai, A 2 , and A 3 denote the eigenvalues whose real 
parts are greater than, equal to, and less than 0 , respectively. 

(Case 1) A > 0. One obtains that the critical coupling 
strength is ref = 27 /A, by setting Re[/z]=0. When re < ref 
[where the linear theory predicts neutral stability. Fig. 0a)], 
r(t) decays exponentially with Re[/r]=Are/2 — 7 =Re[Ai]. 
Therefore, it is remarkable that the damping rate is actually 
determined by the value of the disappeared ghost eigenvalue 
Ai in Fig. 0a), exactly like what happens in the LD mecha¬ 
nism. 

(Case 2) A = 0. The incoherent state is neutrally stable 
regardless of re [Fig. 0c)]. According to Eq. (fTTt . r(t) al¬ 
ways decays exponentially with Re[/i] = — 7 =Re[Ai], once 
again like if the ghost eigenvalue Ai in Fig. 0c) would still 
be present. 




FIG. 3: (Color online) Synchronization transition in Eq. 0 in the 
presence of conformist and contrarian oscillators, (a) r vs A (see 
text for definitions). Letters A, B. C denotes the three conditions that 
will be analyzed in the next Figure. The inset reports r vs p, i.e. the 
proportion of conformist oscillators in the ensemble, which (for any 
given 7 ) is entirely controlled by A. (b) Given 7 , A c comes out to 
be inversely proportional to re. 


(Case 3) A < 0. From Eq. (ITTI) . Re[/r] = Are/2 — 7 , 
which is always negative. Therefore, r(t) does decay expo¬ 
nentially for any value of re, including in the neutrally stable 
regime predicted by linear theory. Remarkably, the decaying 
rate Re[^t] = Are/2 — 7 is here (and must be) different from 
Re[Aa] = Are/2 + 7 = A( K _ because the real part 
of A 3 is positive when re < re]]. Instead, the decaying rate 
Re[/tt]=Re[Ai], i.e., the ghost eigenvalue Ai actually controls 
the decaying of r(t), as shown in Fig. 0d)! 

Finally, we compare our analytical predictions with direct 
numerical simulations of Eq. (0. To that purpose, we set 
N = 10, 000, and adopt a fourth-order Runge-Kutta method 
with integration step At = 0.01. We first verify that the nu¬ 
merics match all theoretical predictions, as shown in Fig. 0 
and confirm an exponential relaxation of R(t ) in the regimes 
of neutral stability, for all the three cases. Next, starting from 
random initial conditions for the phases, we gradually step up 
(forward transition) and down (backward transition) A from 
(and to) a negative initial value (A = —0.2), while keeping 
7 = 0.5 and re = 4.5 constant, this way piecemeal modulat¬ 
ing the proportion of conformist oscillators in the ensemble. 
Results are shown in Fig. 0a), where it is seen that a continu¬ 
ous (fully reversible) synchronization transition occurs at the 
critical point A c = 27 /re. In Fig. 0b), A c vs re is plotted, and 
numerical results perfectly verify the prediction of an inverse 
proportion in the critical point. On its turn, this means that, 
as long as the number of the conformists prevails over that of 
the contrarians, synchronization will always occur for a large 
enough coupling strength. 

Going further, one can even unveil the mechanisms and 
processes taking place during the path to synchronization. 
Namely, Fig. 0 reports the instantaneous distributions of 
phases and frequencies that characterize the coherent states. 
As A is below the critical value, synchronization cannot be 
achieved. When A exceeds A c , a phase-locking cluster of 
conformists (07 > 0) first appears [as shown in Fig. 0a)], 
without an associated cluster of synchronized contrarians. As 
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A increases, more and more conformists join the synchro¬ 
nized cluster. Only at an intermediate moment, the contrar¬ 
ians start forming a synchronized cluster, as shown in Fig. 
|4jb). This latter fact is actually striking, as it demonstrates 
that in such ensembles the synchronization of contrarians can 
only be achieved after the synchronized cluster of conformists 
is large enough, as if the former process would be actually in¬ 
duced by the latter. By further increasing A, the cluster of 
conformists becomes larger and larger by recruiting more and 
more drifting oscillators, which leads to larger order parame¬ 
ter, as shown in Fig. |4] c). The two clusters (that of conformist 
and that of contrarians) rotate with the same frequency along 
the unit circle. During the rotation, they are relatively static 
with each other, as shown in the inset of Fig. fflc2), i.e. the 
two peaks in the phase distribution keep a constant difference 
1.27T (or 0.87t). As a consequence, the system enters (after 
synchronization) into a traveling wave state, which however 
(and generally) is not the 7r state, i.e. that state where the 
phase difference between the two clusters is always 7 r m. 

The occurrence of traveling wave states can be heuristi- 
cally understood as follows. As pointed out in Ref. lEoll . the 
phase-locking condition for a pair of oscillators in system © 
is A ojij = \ijj — Wi\/(\wj\ + |atj|) < nr, i.e. a pair of oscilla¬ 
tors is forbidden to synchronize with each other (for any given 
k ) if such a condition does not hold. Now, if a pair of oscilla¬ 
tors have natural frequencies with different signs, A u>ij = 1 . 
In the incoherent state, nr is typically much less than 1, and 
therefore those oscillator pairs with Aui, j = 1 will most likely 
violate the above condition. At variance, oscillator pairs with 
the same frequency sign are easier to synchronize, especially 
when their natural frequencies are close enough. This explains 
why clusters always form among conformists or contrarians, 
as observed in Fig. [4] In the traveling wave state, both co¬ 
herent clusters of conformists and contrarians rotate with the 
same instantaneous frequency that is greater than 0. Accord¬ 
ing to Eq. ©, the conformists turn to approach the mean-field 
phase f due to uj, > 0. Nevertheless, in order to make the in¬ 
stantaneous frequency 0i > 0 of contrarians, sm((j) — 9) must 
be less than 0, i.e. their phases will always rotate against the 
mean-field phase </>, and one will observe a phase difference 
between the two clusters (conformists and contrarians) which 
is always greater than tt/2. 

In summary, we investigated the synchronized dynamics 
of an ensemble of phase oscillators when the coupling to the 
mean field is frequency-weighted and two populations of os¬ 
cillators (conformists and contrarians) can be identified. Ana¬ 
lytically, we derived the critical coupling strength for synchro¬ 
nization and determined the stability of the incoherent state. 
In all regimes where the linear theory predicts neutral stability, 
the order parameter decays exponentially, in analogy with the 
Landau damping effect in plasma physics. Extensive numer¬ 
ical simulations fully support the theoretical predictions, and 
show that a continuous synchronization transition occurs by 
changing the ratio of conformists and contrarians. Traveling 
wave states are generically observed after synchronization. 
Together with providing generic evidence of Landau damp¬ 
ing effects in the synchronization of coupled oscillators, our 


results are of significance in that they contribute to shed light 
on the mechanisms at the basis of some phenomena beheld 
in social sciences. In many polarized social systems (such as 
two parties in the congress of a country, or fans of two derby 
football clubs) contrarians always confront with conformists. 
In these situations, a phenomenon is frequently observed, i.e. 
the weaker side becomes more united as the stronger sides be¬ 
comes more powerful, while the two sides never compromise 
with each other. Our results may then enlighten the reasons 
for the occurrence of such circumstances. 




FIG. 4: (Color online) Characterization of the coherent states. 

k = 4.5 and 7 = 0.5. Rows (a)-(c) correspond to A = 0.226, 
0.26, and 0.5, respectively, i.e. to the points A, B, and C of Fig. 
[3ja). Column 1 plots the order parameter in the complex plane after 
a transient stage, while columns 2-3 correspond to the snapshots of 
the distributions of the instantaneous phases and the frequencies at 
t = 2, 500. In (a), the conformists slightly prevail over the contrar¬ 
ians, and only a small part of conformists form a coherent cluster, 
rotating at a certain frequency along the unit circle. In (b), the coher¬ 
ent cluster of conformists continuously expands, while the contrari¬ 
ans begin forming a coherent cluster. Both clusters rotate at the same 
frequency along the unit circle, with a constant phase difference be¬ 
tween them. In (c) more conformists are present in the system. The 
cluster of the conformists further enlarges, leading to the increase of 
the order parameter. The inset in (c2) shows the phase distribution at 
that moment. The phase difference between two peaks is 1 . 27 T. 
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